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ANALYSIS OF SUPERSONIC CONICAL FLOWS 


By C. W. Chiang* and Richard D. Wagner, Jr. 

Langley Research Center 

SUMMARY 

An analytical technique is described for conical flow-field predictions for a particu- 
lar class of flows. That is, problems of accelerating and decelerating supersonic cross 
flows are considered; the accelerating cross flow is produced when the free stream is 
normal to the conical axis of symmetry (a fin or stabilizer), and the decelerating cross 
flow is produced when the free stream is in the plane of symmetry (delta wing). For the 
first problem the solution is complete; for the latter only the solution for the supersonic 
cross-flow region is presented. 

The fundamental differential equations are transformed into a dimensionless conical 
coordinate system. Since the partial differential equations governing a conical flow field 
are of hyperbolic type in regions of supersonic cross flow, the method of characteristics 
can be applied. The characteristic equations and compatibility equations are derived. 
Numerical computations are performed, starting from the leading edge and proceeding 
to the central portion of the wing. 


INTRODUCTION 

Linearized theory has been widely used to predict aerodynamic data for slender 
bodies at small angles of attack. For high Mach numbers, especially in the hypersonic 
range, nonlinear effects become important and the results obtained by the use of linearized 
theory are of little value. In fact, at high Mach numbers even the second-order theory is 
not adequate, and meaningful results can be obtained only from analyses including non- 
linear effects. The present work is one phase of a program initiated to develop analytical 
techniques for flow-field predictions of a particular class of flows, taking into account 
nonlinear effects. Only conical flows will be considered. 

Conical flow can be obtained over a conical body; that is, a body whose surface is 
generated by straight lines through a common apex. In conical flow, the flow properties 
along any ray directed from the vertex remain uniform, a feature which allows a trans- 
formation from a three-dimensional Cartesian coordinate system into a two-dimensional 
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conical coordinate system. Thus computations are greatly simplified. Many component 
parts of supersonic vehicles are conical; a flat delta wing is a typical example. 

Numerous examples of the application of linearized theory to supersonic conical 
flows exist in the literature (refs. 1 to 3). The system of rotational conical-flow equa- 
tions has been considered by several investigators (refs. 4 to 15). These analyses show 
that the character of the transformed equations is determined by the cross-flow Mach 
number (i.e., the component of the local Mach number normal to the ray from the apex). 

If the cross flow is subsonic, sonic, or supersonic, the equations are elliptic, parabolic, 
or hyperbolic, respectively. In the present paper, only flows for which the cross flow is 
supersonic are considered. This restriction allows the use of the method of character- 
istics for numerical computations. 

Problems of accelerating and decelerating supersonic cross flows are considered. 
As shown in figure 1, an accelerating cross flow is produced when the free stream is nor- 
mal to the conical axis of symmetry (a fin or stabilizer), while a decelerating cross flow 
is produced when the free stream is parallel to the plane of symmetry of the wing as 
occurs with a delta wing. For the latter problem, solutions for the supersonic cross-flow 
region prevail until the cross flow approaches sonic (though in many cases nearly the 
entire wing is covered); for accelerating cross flows, a complete solution is obtained for 
a symmetrical body. 


SYMBOLS 

a dimensionless sonic velocity (actual sonic velocity divided by Voo) 

aj,a 2 ,ag defined by equations (17) 
bi,b 2 ,b 3 ,b 4 defined by equations (19) 

Cp surface pressure coefficient 

c v specific heat at constant volume 

h enthalpy 

M Mach number 

n unit normal to the surface 

p pressure 

2 


I 


entropy; also shock point 


dimensionless entropy, 


S - S r 


y(y - i)c v 


temperature 


U = — — - 

Voo 

u,v,w 


dimensionless velocity components along X-, Y-, and Z-axis, respectively 
(actual velocity components divided by Voo) 

velocity 


X,Y,Z Cartesian coordinate axes (fig. 1) 


x,y,z 


a 


coordinates along X-, Y-, and Z-axis, respectively 


angle of attack 


a = tan-1 (*»**) 

1 \cos X/ 

0 = cos"l(sin x cos a) 


complement of angle between the surface normal and the stream direction 


ratio of specific heats 


incremental 


gradient vector operator 


flow deflection angle 


e = cj) 1 + a 1 

r = jH— w 
ryu - v 


3 


defined in figure 3 
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£s,max ^-coordinate of upper surface of airfoil at geometrical plane of symmetry 

p density 

<j p i = tan“l( ^ | 

1 (cos x - V sin x) 


(p 2 = tan“l^sin y tan 



X angle between Y-axis and leading edge of airfoil (see fig. Bl) 

to dimensionless vorticity vector 

x, • cross multiply, dot multiply 


Subscripts: 

L leading- edge condition 

N normal to the leading edge 

n normal to the shock 

0 reference condition 

s at the surface 

t tangential 

a behind the shock 

1 first given condition 

00 free stream 


4 



characteristic directions 


+>- 

x,y,z,77,| partial differentiation with respect to x, y, z, 77, and £ 
Accent: 

~ unit vector 


FUNDAMENTAL EQUATIONS 

A rectangular Cartesian coordinate system with origin at the vertex of a delta wing 
is shown in figure 1. The X-axis is directed from the origin along the axis of symmetry 
of the lower surface of the wing. The upper surface of the wing is fixed by the profile 
normal to the X-axis; lenticular and parabolic profiles are used in the numerical calcula- 
tions to be presented. Only attached leading-edge shocks are considered herein, and 
hence the upper and lower wing surfaces are independent. 

Flow Equations 

For steady flow of an ideal gas with a uniform free stream, the continuity, momen- 
tum, and energy equations may be written in the form (ref. 16) 

V • (pU) = 0 (1) 

w x U = VS (2) 

V 1 

v 00 

U • vs = 0 (3) 


where U is a dimensionless velocity vector, actual velocity divided by free-stream 
velocity Voo> 


Flow Equations in Cartesian Coordinates 

The momentum and energy equations (2) and (3) may be written in Cartesian coordi- 
nates in the form 


w(u z ■ 

■ w x) 

- v(v X - Uy) 

= a 2 s x 

(4a) 

u(v x - 

' u y) ' 

- w(Wy - V Z ) 

= a 2 s y 

(4b) 

v(w y ■ 

-V Z ) 

- u(u z - w x ) 

= a 2 s z 

(4c) 
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US X + VSy + ws z = 0 


(4d) 


where 

u,v,w 

s 


dimensionless velocity components along X-, Y-, and Z-axis, actual velocity 
components divided by Voo 


dimensionless entropy, 


S - S r 


In 


P/Po 


y(y - l)Cy y(y - 1) pY/p^y 

S Q reference entropy 

a dimensionless sonic velocity, actual sonic velocity divided by 

Flow Equations in Conical Coordinates 

In dimensionless conical coordinates the position of a straight line directed from 
the apex of the conical body is determined by the coordinates ?j = y/x and £ = z/x. The 
partial-derivative operators in the X,Y,Z axis system may be written in conical coor- 
dinates as follows: 




JL « - 14 , _L + 4 JL 

8x xl 077 s 0£ 


_ 0 _ _ 1 _ 0 _ 
8y _ X 9?7 


_ 0 _ 1 _ 0 _ 
9z “ X 


(5) 




By using the partial-derivative operators of equations (5), equations (4) may be 
written in conical coordinates (ref. 4) : 


= -a. 2 (r}S v + 

(6a) 

= -a^s^ 

(6b) 

) = -a2s ? 

(6c) 

= 0 

(6d) 


Only the last three of equations (6) are independent, and equation (6a) is not used. 
As shown in appendix A, the continuity equation (1) may be written in conical coordinates 
in the form 
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( 7 ) 


2a 2 ^u r? + -v v - w^ - (t?u - v)^ 2 - (£u - w)U^ 2 = 0 
Boundary Conditions 

The boundary conditions are given at the body surface and at the shock (details are 
shown in appendix B). The boundary condition at the body is that the normal velocity 
components vanish. Therefore, 


f dA _ gu - w 
drij s r?u - v 


( 8 ) 


The boundary conditions at the shock can be found by applying the Rankine-Hugoniot 
relations across the shock. To facilitate the application of the shock conditions, a ficti- 
tious plane shock which is tangent to the body shock is introduced. The flow properties 
across this shock will match those behind the body shock at the point of tangency. In a 
plane which is normal to the fictitious leading edge, and hence the fictitious shock plane, 
the angle e between the fictitious shock plane and the plane of the undisturbed flow may 
be determined for known values of the free-stream Mach number Moo, the angle of attack 
a, the flow deflection angle 6, and the sweep angle y, from the following equation: 


1 + 


r - 1 


M 0 


sin 2 ^)tan 3 e - (Moo 2 sin 2 /3 - l) cotta 


tan 2 e 


+ 1 + 


y + 1 


M 0 


tan^tan e + cot^ + 6j = 0 


( 9 ) 


Once the angle e is determined, the velocity components behind the shock and the 
entropy change across the shock can be found. 

The velocity components behind the shock for the problem of decelerated cross flow 
(wing problem) are 


u = cos /3 sin y + 


cos y cos 5 sin /3 cos e 


cos (t 


- “l " 6 ) 


v = -cos )3 cos y + 


sin y cos 6 sin /3 cos e 


cos^e 


- “i - 6 ) 




w = 


sin 6 sin /3 cos e 
cos^e - tXj - 6' 


( 10 ) 


Similarly, for the problem of accelerated cross flow (fin problem), the velocity 
components behind the shock are in the form 
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u = -cos /3 sin x + 


v = cos /3 cos x + 


cos x cos 6 sin ft cos e 
cos^e - CKj - 6) 

sin x cos 6 sin ft cos e 
cos(e - - 6) 


_ sin 5 sin ft cos e 


cos^e - aj - 6j 

The entropy change across the shock is calculated from 




Vs = s - s c 


y(y - 1) 


In 


^- (Moo^sin^j3 sin^e - — ~ ^ 

y + 1\ y + 1/ 


In 


(y + ljMoo^sin^ sin^e 


2 + (y - ljM^sin^/S sin^e 


( 11 ) 


( 12 ) 


CHARACTERISTIC AND COMPATIBILITY EQUATIONS 


Characteristic and compatibility equations are derived in detail in appendix C. In 
addition to boundary conditions, they constitute the necessary requirements for actual 
computations. 

Along a stream surface a characteristic is given by 

d| = ju_-_w (13 

d-rj tju - v 

For this characteristic there is one compatibility relation due to constancy of 
entropy along stream surfaces: 


ds = 0 


(14) 


The second compatibility equation along stream surfaces may be written in the 


form 


du + r] dv + ij dw = 0 (15) 


The remaining two characteristics, one above and one below the stream surface, 
have slopes given by 


/dA a l ± fl ~ a 2 a 3 


IdTJi 


a 2 


(16) 
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where 


a-i = a 2 £rj - (tju - v)(£u - w) 
a 2 = a 2 (l + tj 2 ) - (tju - v) 2 ) 
a 3 = a 2 (l + £ 2 ) - (£u - w) 2 ^ 

The compatibility equation corresponding to the two characteristics is 

bj du + b2 dv + b 3 dw + b4 ds = 0 


where 



■\ 


(|u - w) - (77U - v]^- 


b 2 = a 2|( 

b 3 = -a 2 £(£v - rjw) - £u|ai - a 2 + 


j - a 2 ?7(^v - rjw) - ^u^ - a 2 


a 3 a 2 


(|u - w) - (tju - v)|| 


b^ = -a 2 


f dA 

a i " a 2 -r) 

v 1 <W 


Q 

a2 d^ 


( 17 ) 


(18) 


(19) 


METHOD OF COMPUTATION 

The application of the method of characteristics to compute conical flow fields con- 
sists of the repeated use of three basic computing processes (one each for the shock 
point, the body point, and the general field point) proceeding inward from the leading edge 
with given initial conditions. With reference to figure 2, the initial conditions are estab- 
lished by approximating the leading portion of the wing by a swept wedge, tangent to the 
wing at tjj. The oblique-shock solution, equations (12) and (10) or (11), then gives the 
flow conditions at the initial station ?7j. 


Shock Point 

To determine the ij- coordinate of the shock at tj = ?jj + AT7, as a first approxima- 
tion the shock at Tjj is extended to Tjj + Atj to establish the point S. The shock rela- 
tions give the first approximation to the flow properties at S, and the characteristic SAj 
establishes the point Aj. The flow at point Aj is known by interpolation, and the com- 
patibility relation, equation (18), must be satisfied for the characteristic SAj as well as 
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the shock conditions at S. An iteration on the shock angle at S, with averaged coefficients 
in the compatibility equations, allows determination of the correct shock angle and the 
correct locations of S and Aj, and thence the flow properties at point S. 

Body Point 

The point A p on the characteristic A p P p intersecting the body at 77 = ?7j + A77 can 
be found by interpolation. The compatibility relation for A p P p , the two compatibility 
relations for the streamline characteristics (eqs. ( 14 ) and ( 15 )), and the boundary condi- 
tion (eq. (8)) can then be solved for the flow conditions at P p . 

General Field Point 

The present calculations were performed by computing conditions from the shock to 
the body. After the point Aj is established, the characteristic AjPj locates the first field 
point Pj. The streamline P^B^ and the characteristic CjPi through P^ can be determined 
by interpolation of the given, or known, conditions on 77 j. The compatibility equations for 
the characteristics AjP^ and CjPi and the two compatibility equations for the streamline 
PjBi can then be solved for the flow conditions at point Pi. This process is repeated for 
the point P2, and so on, until a point Pj is reached with the characteristic PfCj intersecting 
the initial data line r\ - 'Qj (or the current known data line) below the body surface. The 
body point is then calculated; the flow conditions are then completely specified along the 
line 77 = 775 + A77, from the shock to the body. This new set of data can be used to establish 
the conditions at the next station, v - Vi + 2A?7. The computations thus proceed over the 
wing, or fin, until the sonic condition is obtained, or the trailing edge is reached. 

The method as presented with iterations and quadratic interpolations for the flow 
properties is accurate to second order (i.e., the truncation error is the order of the mesh 
size cubed). This computational method was programed in FORTRAN language for com- 
puting on an electronic digital computer. 

APPLICATIONS 

Accelerating Cross Flows (Fin Problems) 

Calculated pressure distributions for fins with circular-arc profiles are shown in 
figure 3 . The results are given in terms of the ratio of the local surface pressure coef- 
ficient C p to the maximum pressure coefficient C p? L, which occurs at the leading edge. 
The fin geometry is defined by the sweep angle x and the initial flow deflection 6jsf in 
a plane normal to the leading edge. Calculations for Moo = 6, 10 , and 30 were performed 
for fins with y = 30 ° and = 5 ° and 20 °. The results are presented up to the geomet- 
rical plane of symmetry of the fin, 6 = 180 ° (0 is defined in fig. 3 ), although actual cal- 
culations were carried beyond that. Calculations can be carried close to the trailing edge 
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with much finer meshes and longer computing time. Included in figure 3 are calculations 
based on a shock-expansion method for conical flows which was obtained in reference 17. 
The agreement of the two methods is exceptional, considering the simplicity of the shock- 
expansion approach. Because earlier studies have shown that the two-dimensional shock- 
expansion method works best at high values of y, a cursory examination of the effect of 
y was made. In figure 4 calculations are shown for a fin with a circular-arc profile, 
ot — 0°, = 30°, X = 60°, and = 30 at y = 1.667, 1.400, and 1.100. The shock- 

expansion method works best for y = 1.400 and the poorest agreement is for y = 1.100. 
These results are consistent with the y-effect observed in the prediction of two- 
dimensional flows by two-dimensional shock expansion. 

For two-dimensional and axisymmetric flows, the application of the shock-expansion 
method is greatly simplified because of a closed-form relation between the surface pres- 
sure and the local surface inclination. For conical flows such a relation is not known and 
solutions can be obtained only by quadrature; more easily obtained approximate solutions 
are desirable. Newtonian theory is often used by the practicing engineer in the analysis 
of complex three-dimensional configurations. There are too few exact solutions to enable 
an assessment of the applicability of Newtonian theory to three-dimensional hypersonic 
flows. Indeed, there has been no assurance that the usefulness of Newtonian theory for 
two-dimensional and axisymmetric flows should prevail for three-dimensional flows. In 
figure 5, results obtained from Newtonian theory have been included for comparison with 
results obtained by the method of characteristics, the shock-expansion method, and a strip 
theory. The last of these is a two-dimensional approximation using the circular profile 
and two-dimensional shock-expansion theory. The Newtonian theory used is a logical 
extension of the generalized Newtonian theory (ref. 18) and the pressure ratio is given by 

Cp _ sin 2 r 
C P,L sin 2 r L 

where r is the complement of the angle between the surface normal and the stream 
direction, and Cp^ is the pressure coefficient at the leading edge, where r = T-^, 
obtained from the oblique-shock relations for the tangent swept wedge. The ratio of the 
surface pressure to leading-edge surface pressure was calculated at y = 1.667 and 1.400 
for a fin with Sjq- = 20°, x = 30°, and = 30. Strip theory always underestimates the 
pressure and the Newtonian theory gives poor results in the low-pressure regions (a gen- 
eral failing of Newtonian theory; see ref. 19). 

In figure 6 are shown comparisons at Moo = 6.85 between the method of character- 
istics, shock-expansion method, and unpublished experimental data obtained by M. H. 
Bertram at the NASA Langley Research Center. In these tests the ratio of wall tempera- 
ture to total temperature was in the range of 0.6 to 0.8, the free-stream Reynolds number 
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was 0.31 x 106 per inch (0.12 x 10® per cm), the root-chord length was 4.8 inches 
(12.2 cm), and the semispan was 2.4 inches (6.1 cm). The local chord was 2.38 inches 
(6.05 cm) at the section in which the pressure orifices were located. The profile of the 
surface consists of two parabolas in the YZ -plane, tangent at r) = \|3/2; namely, 

9 

ij = q tan xiv cot X - 0.5) 

where q = 0. 1 for ^3/2 < r/ = ^3 and q = 0.05 for 0 = 77 < ^3/2. The maximum 
thickness occurs at 17 = ^3/2. In figure 6 , x = 30°, and ~ 11 ° at a = 0°. The two 
theoretical solutions agree well, and the difference between the experimental data and 
theory can be attributed to the viscous boundary-layer displacement effects (ref. 20). At 
a = 10 ° the theory and experimental results are quite close. 

Decelerating Cross Flows (Wing Problems) 

In the case of a delta wing where the free-stream velocity vector lies in the plane 
of symmetry, the cross flow is supersonic at the leading edge and decelerates to sub- 
sonic and ultimately to zero at the plane of symmetry. The present calculations can be 
used only in the region of supersonic cross flow, but this flow region can constitute a 
major portion of the flow field. (For example, for a circular-arc wing with Moo = 30, 

6 j^ = 20°, and x = 30°, about 90 percent of the wing flow field can be computed.) 

Calculations for a circular-arc wing are shown in figure 7. The wing parameters 
are = 9° 18' and x = 50°. The shock shape and the pressure ratio of both the shock 
and the wing surface are given for Moo = 4.0 and a = 0°, 5°, and 10°. For this case 
the method of characteristics can be used to compute a major portion of the flow field, 
until the cross flow approaches sonic as shown in figure 7(a). As the angle of attack 
increases, the portion of supersonic cross flow decreases. A line of constant cross-flow 
Mach number of 1.08 for the case of a = 0° is also shown in figure 7(a) to give some 
idea of the Mach cone which separates the supersonic and subsonic cross flows. In ref- 
erence 21 calculations based on the method of lines, giving a solution for the entire flow 
field, are presented for the same conditions; these are shown in figure 7(b). The agree- 
ment of the two methods is excellent. The pressure ratio calculated on the basis of 
Newtonian theory is also shown; it has lower values in the supersonic region than the 
other two methods. 

An additional comparison of the method of reference 21 with the method of charac- 
teristics is given in figure 8 , which shows the surface pressure coefficient on a circular- 
arc wing with x = 50°, = 20°, a = 10°, and M^ = 8.1. Also included in figure 8 is 

the shock- expansion method. All the methods show excellent agreement. The shock 
shapes predicted by the method of characteristics and the method of lines also are 
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compared and they agree closely. In reference 21, the solutions by the present method 
and the method of lines are compared and verified with experimental data. 

CONCLUDING REMARKS 

An analytical technique using the method of characteristics has been described for 
flow predictions of accelerated and decelerated supersonic cross flows over a conical 
body. Results obtained by this method agree well with experimental results and other 
analytical results. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 28, 1970. 
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APPENDIX A 


CONTINUITY EQUATION IN CONICAL COORDINATES 


The entropy function can be written in the form 


= e^r-!) 3 =f( s ) 


P y Poo 7 


and enthalpy may be written in the form 


h = 


= y P - 1 a 2 

y - 1 P v „ 2 y - 1 


Combining equations (Al) and (A2) gives 


where 


p = k 1 (a 2 )v~l|f(s)J r- 1 

1 


k l = 


2W-1 


When equation (A3) is substituted into equation (1) in the body of the paper, the 
tinuity equation becomes 

‘ JL -11 

V • (pU) = kqV • [a y_1 f(s) y - 1 UJ 

— 1 

= kjf(s) y “^V • _a y ~^UJ - ■ - j f(s) y ”^|a y_A Uj • Vf(s) 

= 0 

Since U • VS = 0, then 

U • v[f(s)] = y(y - l)f(s)U • Vs = 0 


A\Av] ■ 


Thus 


but 


2 ' 
a^U 


= 0 


a 2 = (r-l(h 0 -^) 


(Al) 

(A2) 

(A3) 

con- 


(A4) 

(A5) 
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APPENDIX A — Concluded 


and combining equations (A4) and (A5) gives 

2a 2 V • U - U • VU 2 = 0 

Equation (A6), the continuity equation, may be written in conical coordinates as 

2a 2 (v u n + " v 77 " w ^) " ( r 7 u “ v ) u tj 2 " “ W )U| 2 = 0 


(A6) 


(A7) 
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APPENDIX B 


BOUNDARY CONDITIONS 


The boundary condition at the body is that the normal velocity components vanish. 
If the equation at the conical surface is g (£, 77 ) = 0, then the boundary condition at the 
body is 


or 



d|\ _ £u - w 
*7/ s ~ m - v 


(Bl) 


The boundary conditions at the shock can be found by applying the Rankine-Hugoniot 
relations across the shock. The equations presented herein are extensions of those given 
in reference 4 for the conditions at the leading edge of a flat delta wing; they are general- 
ized here to apply to any location along the shock contour. To help in understanding the 
relations in front of and behind the shock plane, a geometrical layout is constructed as 
shown in figures Bl and B2. To facilitate the application of the shock conditions, a 



Figure Bl.- Geometrical layout. 
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APPENDIX B — Continued 


fictitious plane shock which is tangent to the body shock is introduced (plane OFM in 
fig. Bl). The flow properties across this shock will match those across the body shock. 
The equation of the fictitious shock plane is (by geometry) 

tan <f> i = (B2) 

cos x - V sin x 

where 0j is the slope of the fictitious shock plane measured in the plane which is nor- 
« mal to the fictitious leading edge (in plane FHG). 

The angle between the direction of the undisturbed flow and the fictitious leading 
% edge, /3, and the angles and cp 2 may be obtained through analysis of the geometry: 

/3 = cos“l(sin x cos a) (B3) 




(b) 

Figure B2.- Velocity components with rotating axes. 
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APPENDIX B - Continued 


«, = tan-ipS^L] (B4) 

1 \COS X/ 

tan 02 = sin x tan 0 (B5) 

where 0 2 is the slope of the shock plane in the plane normal to the X-axis. 

The plane FHN is redrawn in detail in figure B2(a). The free-stream velocity com- 
ponent in the direction of HN is V*, sin jS, and the normal and tangential free-stream 
velocity components are designated as V n>0 o and Vt ?0 o, respectively. The resultant, 
normal, and tangential velocity components behind the shock are denoted by V a , V n , and 
Vj-, respectively. 

At the leading edge of the fictitious -plane shock, the shock wave must satisfy the 
relation (ref. 22) 


V„ ,« • V„ = (a*) 2 - V t 2 

where a* is critical sonic velocity. 

From figure B2 it is easily seen that 

V t, OO — Voo sin /3 cos e 
V n ,°o = Voo sin /3 sin e 


(B6) 


(B7) 

(B8) 


where e = 0j + represents the angle between the fictitious shock plane and the plane 
of undisturbed flow. 

Since 

Vt = V t ,oo 


^n = v v 


tan^e 


a l) = V 00 sin /3 cos e tan^e - 6 - a-j^ 


(B9) 


where 6 is the angle between the cross-flow velocity and the XY-plane measured nor- 
mal to the fictitious leading edge. 


The energy equation may be written 


h 



+ 



2 


or 


aoo 2 + ^(Voo sin = (a*) 2 + *-^(a*) 2 = a*) 2 
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APPENDIX B - Continued 


(a*) 2 = — a J + ( V - sin Z 3 ) 2 (B : 1°) 

y + 1 y + 1 7 

where V* and a* are critical velocity and critical sonic velocity. 

Combining equations (B7), (B6), and (BIO) yields 


Voo^sin 2 /3 sin e cos e tan/e - 6 - q'i') = — aoo 2 + - - \ Voo 2 sin2/3 - y - \ Voo 2 sin^/3 cos^e 

V V y + 1 y+1 y+1 

If tan(e - 6 - a is expanded, Vco/aoo is replaced by M*, and both sides of the equa- 
tion are divided by tan (6 + a-^, the result is 


1 + y ~ Moo 2 sin2^Jtan3e - ^Moo 2 sin2/3 - ljcot^5 + a^tan^e 
1 + ^ M^sin^tan e + cot ^5 + a-jJ = 0 (Bll) 


For known values of M^, a, 6, and x> equation (Bll) determines the angle e; 
the velocity components behind the shock plane and the entropy change across the shock 
can then be calculated. 

To obtain velocity components behind the shock, the reader is again referred to 
figure B2. The coordinate axes j’,k' and j",k" are parallel and normal to HG and HJ, 
respectively. The coordinate axes i,j are parallel to the X-axis and Y-axis; the axis i' 
is in the direction of the fictitious leading edge, and the axis j' is normal to the axis i' 
and is in the direction of HG, which coincides with the XY-plane (see fig. Bl). 

By coordinate transformation, the unit vector in the direction of the velocity V CT 
behind the shock is 

j" = i cos x cos 6 + j sin \ cos 6 + k sin 6 (B12) 


V sin g cos e (B13) 

cos(e - 6 - a-d 

Equations (B12) and (B13) then define the magnitude and the direction of the velocity 
V a behind the shock. The total velocity is simply the vector sum of this velocity vector 
and the free-stream velocity component tangential to the leading edge, 

V ro cos fi(i sin x - j cos x). 

The velocity components behind the shock for the problem of decelerated cross flow 
(wing problem) are then 
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APPENDIX B - Concluded 


u = cos 0 sin x + 


cos x cos 5 sin 0 cos e 
cos (e - 6 - 


v = -cos 0 cos x + 


sin x cos 6 sin 0 cos e 


cos 


( € * 6 “ “l) 




(B14) 


_ sin 5 sin 0 cos e 


cos 


(e - 6 - a x ) 


Similarly, for the problem of accelerated cross flow (fin problem), the velocity 
components behind the shock are 


u = -cos 0 sin x + 


cos x cos 5 sin 0 cos e 


v = cos 0 cos x + 


cos^e - S - a-jj 
sin x cos 6 sin 0 cos e 


cos 


( e - 6 " «l) 


(B15) 


w = 


sin 6 sin £ cos e 


cos 


(e - 6 - aj) 




The entropy change across the shock is calculated from (ref. 22) as 


As = s - Soo = 


y(y - l) 


In 


y + 1 


^—(Moo^sin20 sin2e 


y - l 

y + 1 


y - 1 


In 


(y + l)M oo ^sin^0 sin2g 


2 + (y - l)M oo 2 sin 2 0 sin 2 e 


(B16) 
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APPENDIX C 


CHARACTERISTIC AND COMPATIBILITY EQUATIONS 

Any conical surface may be expressed as a function of the two coordinates g 
and 77: 


g(g,77) = Constant 
dg = dg + g^ dr] = 0 
The unit normal to the surface is 


n 


= _Zi_ 

Vg 


Maslen (ref. 15) shows that the cross -flow velocity component — that is, the component 
normal to a ray from the apex — governs the type of partial differential equations in con- 
ical flow. For supersonic cross flow he shows that the characteristic surfaces are 


U • n = a and U • n = 0 
where U • n can be expressed as 


U • n = 




t dg \ 2 /dg\ 2 n 

+ y +i 


(Cl) 


For the case where U • n = 0 the characteristic surface is given by 


dg _ gu - w 
d77 77U - v 

which is a surface consisting of streamlines. 

From equation (6d) in the body of the paper, the compatibility relation is 


(C2) 


ds = 0 (C3) 

This, of course, expresses the constancy of entropy along streamlines. The character- 
istic (C2) is a double characteristic since multiplying equation (6b) by rju - v and equa- 
tion (6c) by gu - w and adding them gives 

21 


1 



APPENDIX C - Continued 

0?u - v)^Ujj + riv v + £w^ + (£u - w)(u| + ? 7 V| + £w^ = 0 (C4) 

It is found by inspection that equation (C2) defines a characteristic for equation (C4) for 
which the following compatibility relation applies: 


du + 77 dv + £ dw = 0 


(C5) 


The equation U • n = a is quadratic, and the two roots yield the slopes of two 
characteristics, one above and one below the stream surface: 



H * V a i 2 - a l a 3 


a 2 


(C6) 


These are the projections of the Mach cones in the cross-flow direction. 

The remaining compatibility relations are found as follows. Equation ( 6 d) and 
ds = Sjj dr) + s^ d| yield the partial derivatives of s in terms of total differentials: 


where 


s t =• 


ds 




d? - ? d77 


S„ = — 

77 d|-?d r? 

? _ |u - w 

7 7U - V 


(C7) 


(C8) 


Multiplying equation (C4) by dr] and then adding and substracting the quantity 
(u| + ?}V£ + £wAd£ yields 


where 


u^ + rjV| + £w^ = Fjj 
u^ + r]V v + Z-w v = F 2 ) 




Fi = 


1 d£ - ? d?? 
?P 


F 2 = 




d£ - ? dTj 
P = du + r) dv + £ dw, 


(C9) 


(CIO) 
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APPENDIX C — Continued 


Substituting equations (C9) into either equation (6b) or equation (6c) yields 


where 


Q = 


V £ - w 7? = Q 

a 2 ds + uP 


(t?u - v)(dij - £ dii) 

The continuity equation (7) may be written in the form 

a 2 £u£ + a 2 7ju 7? + (r/u - v)(|u - w)v| - |a 2 - (?7u - v) 2 j 

- |a 2 - (£u - wj^jw^ + (tju - v)(£u - w) = 0 


or 


where 


Since 


a l v 4 + a 2 v ?7 + a 3 w ^ + a l w ?j = G 


G = -a 2 ^Fi + 77F2) 
aj = a 2 £?7 - (tju - v)(£u - w) 
a 2 = a ^(l + 17 2 ) - (nu - v) 2 
a 3 = a 2 (l + £ 2 ) - (£u - w) 2 

J 


(Cll) 


(C12) 


(C13) 


(C 14) 


V| d| + Vjj drj = dv 
W| d| + w^ dr) = dw 

equations (Cl 5), equation (Cll), and equation (C13) can be written in matrix form as 


(C15) 


where 


Ax = B 



a i 

a 2 

a 3 

a l 




G 


1 

0 

0 

-1 


V- 


Q 

A = 





X = 

V 

B = 



d| 

dt) 

0 

0 




dv 


0 

0 

d| 

di) 


Wrj 


dw 


(C 16) 


(C17) 


23 



APPENDIX C - Concluded 
The determinant of the coefficient matrix A is 

| A | = d£ 2 - 2aj d| drj + a.% d? 7 2 (C18) 

When the determinant vanishes, it gives the slope of the characteristics, equation (C6). 

To find the compatibility equations, the determinant 


a l 

a 2 

a 3 

G 

1 

0 

0 

Q 


d?7 

0 

dv 

0 

0 

d? 

dw 


is equated to zero, or 




dv 


dw 


G . a lQ+ _ a2 Q.. a 2 _.a3_ = 0 


Substituting values of G and Q gives 

bj du + b 2 dv + b 3 dw + ds = 0 


where 


t>i = -a 2 (£v - ? 7 W) - /aj - a 2 ~ju 


b 2 = a 2 


(|u - w) - (rju - v)|i 
drj 


- a 2 rj(£v - ryw) - 7711 /ai - a 2 ^ , 


b 3 = -a 2 £(£v - rjw) - ^ufaj “ a 2 — jj + • 


a 3 a 2 


(?u - w) - (rju - v> 


d| 

dr] 


d| 

drj 


b4 = -a 2 L - a 2 m 


(C19) 


(C20) 


(C21) 


Equation (C20) is the compatibility equation, which is the condition to be fulfilled along 
any characteristic given by equation (C6). 
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Figure 3.- Accelerating cross flows over circular-arc fin. x = 30°; y = 1.400; a = 0°. 




Figure 4- Effect of y for circular-arc fins. Mco= 30; 6fj = 30°; x = 60°; a = 0°. 
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Figure 7.- Shock-wave shape and pressure ratio as a function of percentage of half-span for a circular-arc wing. 

r = 1.40; Moo = 4.0; S N = 9018’; x = 50°. 





— Method of lines (reference 21) 
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